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Abstract 

Recent experimental results for the reaction y + p — > K + + 2 + n from CLAS at Jefferson 
Lab are discussed. It was found that the mass distributions or "line shapes" of the three charge 
combinations £ + 7r~, E°7r° and £~7r + differ significantly. Our results show that the A(1405), as 
the 7 = constituent of the reaction, must be accompanied by an I > component. We discuss 
phenomenological fits to the data to test the possible forms and magnitudes of these amplitudes. 
A two-amplitude I = fit of Breit-Wigner form to the iPn channel alone works quite well. The 
addition of a single 7=1 amplitude works fairly well to model all the line shapes simultaneously. 
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1. Introduction 

The CLAS Collaboration at Jefferson Lab has measured [1|, as part of a large program of 
photoproduction reactions in the GeV energy region, the channel y + p — > K + + £ + ji. The Ztt 
invariant masses range was from very close to threshold at 1328 MeV to well beyond 1520 MeV 
where the A(1520) is prominently produced. In between these limits, the A(1405) is produced, 
as well as the 2(1385) and perhaps other possible states. The CLAS experiment is the first to 
simultaneously map out this mass range in all three charge combinations £ + 7r~, E°7r° and 5T7r + . 
This is interesting because comparison can reveal (1) whether there is substantial contribution 
of / > amplitudes to the reaction mechanism even after the 2(1385) is removed, and (2) how 
strong the channel coupling is between the Ett final state and the competing NK final state. Thus, 
while there have been hints of line shape distortions of the A(1405) in earlier experiments [2-8 1 
the CLAS experiment offers the best look so far at this phenomenon. 

Satisfactory theoretical understanding of the A(1405) continues to be a challenge. For a 
recent review, see Ref 0. For example, within the spectrum of P-wave baryons within the 
relativistic quark model it is a poor quantitative fit liTOl ITT1 . In models of so-called chiral unitary 
state generation among octet baryons and pseudo-scalar mesons, this state plays an important 
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Figure 1: (Color online) Creation of the 
three-body K + £ n final state via an intermedi- 
ate hyperon in the reaction y+ p — > K + +1.+7T. 
In this particular example, a /-channel ex- 
change enables an off-shell kaon to create a 
A(1405) that is sub-threshold for on-shell NK 
reactions. 



role lITSl IT3l . It arises through channel-coupling of the NK system, with threshold at 1432 MeV, 
and the Ym system into which the A(1405) must decay. In the flavor symmetric SU(3) limit, 
chiral unitary models show lfT4ll that two new octets and a singlet of J p — i states are formed 
which, when SU(3) is broken, leads to the presence of two 1 = poles in the hyperon sector 
that reside near the mass of the A(1405). Meanwhile, the 1=1 poles vanish or are very weak. 
Thus, it is of interest to see how well the "two-pole" picture of the A(1405) is supported by 
the data. One prediction is that the line shape of the A(1405) depends on the reaction channel 
through which it is excited, since a given reaction can couple with differing strengths to the two 
predicted poles. In the present case we look using real photoproduction off a proton target. This 
is illustrated in Fig.[T]which shows how an off-shell kaon exchange can create the A(1405) that 
is sub-threshold for a free kaon plus a nucleon. A calculation of this channel was reported in 
Ref. |fT5ll and compared to our results in Ref. (TJ. 

The CLAS measurements are described in detail in Ref. [1|. Briefly, tagged, real, unpolar- 
ized photons [16| with a highest energy of 3.8 GeV impinged upon a 40 cm liquid hydrogen 
target. About 20 billion events with at least two charged tracks in the large-acceptance toroidal 
spectrometer ifTTl were recorded, and in the analysis about 1.41 x 10 6 events with K + Y,n were 
isolated. Charged kaons, protons and pions were identified by time of flight over a four to five 
meter flight path. 

The 2 + was reconstructed via both of its decay modes, a proton with an undetected 7r°, or a 
n + with an undetected neutron. The IT was reconstructed via its decay to if and an undetected 
neutron. Regions of kinematic overlap of the two charged Z's were cut away. The iPn final 
state was reconstructed through its decay via A — > n~p with an undetected n° and photon. The 
undetected particles were identified via missing mass, and a 1C kinematic fit was employed to 
optimize resolution. The exception was the iPn case where two neutral particles were missing, 
so the selection was made using a fixed range of missing masses. 

The two identified backgrounds removed were from ^ + 2°(1385) production and from K*lP 
production. The first of these was measured directly in its dominant decay mode to At: in P- 
wave. Using the 1 1.7% branching fraction to Yfn* we could accurately subtract this contribution; 
this was generally a small fraction of the total data. Coherent interference with S-wave decays 
was in principle canceled due to the full angular coverage of the Y.ti final state in the Y* rest 
frame. K*T. production resulted in a broad background underneath the whole of the ~Zn spectrum, 
with varying kinematic overlap in different bins of W. Studies showed that making increasingly 
drastic rejection cuts of events in the wide K* region made no change in the line shapes of the 
~Ln results. We were unable to find any hint of coherent interference due to this background. 
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Both backgrounds were modeled in a Monte Carlo simulation and subtracted away. Details of 
the analysis procedures, including systematic studies, can be found in Refs. |[T| and Ifl8l . 

The analysis used 100 MeV wide bins in yp center of mass energy W to present the results, 
and it will be seen that the line shapes differ as W increases. 

In this article we present two fits to the "Zn line shape data that are different from the "best" 
overall fit described in Ref. (T). In particular, we first show the fit obtained when allowing two 
7 = Breit-Wigner (BW) amplitudes to fit only the L°7r° data. This final state cannot be reached 
via an I = 1 decay by reason of the isospin-addition coefficients, and we neglect 7 = 2 processes. 
Thus, this fit can serve to test the two-pole picture of the A(1405). Second, we present a fit that 
adds a single 7=1 BW amplitude to the previous two 7 = amplitudes and use this to fit all 
three Y.n states at once. This serves to characterize the isovector piece of the interaction needed to 
represent the experimental data. This amplitude combination leads to a result that is statistically 
less good than the one described in Ref. [1| which used one 7 = amplitude and two 7=1 
amplitudes. On the other hand, it is closer to some theoretical biases about the likely two-pole 
structure of the A(1405). Section[2]outlines the fit procedure, and Section[3]discusses the results. 



2. Model for Mass Distributions 

As described in Ref. [1J, we posit that the E^-part of the three-body K + , S, n final states arises 
via an interaction schematically labeled T^'\ where 7 denotes the isospin of the "Ln system. Let 
the matrix element squared be denoted by 

|f / | 2 = |<7,7 3 -0|f (/) | r p)| 2 . (1) 

The matrix element expressions that give the three measured final state charge combinations are 
then given by the coherent sum of the appropriately Clebsch-Gordon-weighted tj elements as 
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At a fixed yp center of mass energy W the matrix elements are parametrized as relativistic Breit- 
Wigner functions Bj{m), where m is the "En mass, with an overall strength C(W) and a strong 
production phase A<f>i as 

f/(m) = C / (W)e iA0 '7i / (OT). (5) 



The Breit-Wigner amplitude form is 

Bi(m) = 



m 2 R -m 2 - im R T tot (q) 



(6) 



where mg is the centroid of the resonance distribution, T~ is the fixed decay width to a given final 
state, and r tot (g) is the total width to all final states. The available momentum in the decaying 
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hyperon center of mass system is q, and in this frame q R is the available decay momentum at 
m = m R . Angular momentum L is zero in the present case of an S-wave decay of the odd-parity 
A(1405) decaying to a pseudo-scalar meson and an octet baryon. 

The mass-dependent width in the denominator was treated in a way that accounts for channel 
coupling using the Flatte prescription lfl9l . That is, if decay mode "1" is the Y,n final state, and 
decay mode "2" is the NK final state, then part of the mass range is below threshold for mode 
2. Nevertheless, mode 2 has influence on mode 1 both below and above threshold. To preserve 
unitarity and the analytic form of the amplitude, we analytically continuing the center of mass 
momentum q of decay mode 2 to below its threshold. We write 

TtotO") = r 7i i(^i(m)) + r It2 (q2(m)), (7) 

where the decay channels are described by width 

rz,X,) = ^f^r. (8) 

' J m \ q R ) 

Below mode 2 threshold mthresh; the momentum qiim) is nominally zero. However, we continue 
the momentum to imaginary values for m < m t h res h- Furthermore, we introduce a Flatte factor for 
the branching fraction of the decay modes as 

r = rt 2 /r°i- (9) 

This factor is determined by the fits, and indicates the ratio of the channel coupling. 

Finally, the mass distribution is computed as a cross section differential in "En mass m by 
including all the flux and phase space factors for the initial and final states as 

do-ab (hc) 2 a p K +q 2 

~d^= 64.3 PypW 2^ ■ W 

Here p K + is the cm. kaon momentum and p yp is the initial state cm. momentum. We have 
factored out the electromagnetic coupling strength a so that C(W) becomes an effective strong 
coupling in units of VGeV. This expression is integrated over all kaon production angles, a step 
that was necessary to have adequate statistics. Similarly, it is integrated over all 2 decay angles 
in the hyperon rest frame, as mentioned, to ensure cancellation of any interference with P-wave 
decays from the 2(1385) (which is anyway known to be small and was incoherently subtracted). 



3. Results 

First we consider the fit obtained when using two 7 = amplitudes fitted to the data for 2°7r° 
only. This decay mode should be dominated by the 1 = amplitude, assuming 1 = 2 can be 
neglected. Figure [2] shows a sample result in one 100 MeV wide bin of W at 2.3 GeV. The solid 
blue curve shows the total fit, and the solid and dashed black curves show the cross section that 
would result from each of the two amplitudes separately. Also shown is a linear background (not 
part of the fit) from threshold to 1 .6 GeV that was included to account for the imperfect modeling 
of background in the experimental analysis. Data points near the location of the A(1520) position 
have also been suppressed. The fit was done not just to the W bin shown, but to all nine bins at 
once. The full data set is shown in Fig. [3] In each panel the black curves have the same shape, 
differing only in overall weight. The weights, the C/(W), are shown in Fig. |4] 
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Figure 2: (Color online) Two / = ampli- 
tudes fitted to only (solid and dashed black 
curves). The error bars are combined statistical 
and systematic point-to-point uncertainties. The 
solid blue curve shows the fit to the lPn° data 
(blue points). The three vertical dotted lines at 
the left are thresholds for the three decay modes, 
and vertical black dotted lines mark the nomi- 
nal 1 .405 GeV mass and the location of the NK 
threshold m t h re sh- The incoherent background is 
shown as a thin dashed line (blue). 
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It is evident visually that the fit is good; the overall reduced^ 2 of the fit was 0.89. The results 
of the fit are given in Table [T] The two 7 = amplitudes have centroids at 1329 and 1390 MeV, 
respectively, with respective widths that both differ when compared to the 50 MeV RPP value 
for the A(1405) width. This is in part explained by the action of the channel coupling. It is 
seen in the table that the Flatte factor for the narrower amplitude (the only one permitted in the 
fit) is large, at 1.5. This has the effect of distorting the line shape around the NK threshold, as 
seen in the figures, but it also has the effect of making the most probable mass somewhat above 
1405 MeV, even though the nominal centroid is at 1329. The goodness of this fit may be taken 
as positive evidence for the two-pole solution for the 7 = state(s). Figure [4] shows how the 
two components of the fit vary in strength as W increases: evidently the broader one fades out 
quickly while the dominant narrower one remains strong. We have no interpretation to offer for 
this behavior, nor for the relative production phase between the two amplitudes given in the table. 
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Figure 3: (Color 
online) Two 7 = 
amplitudes fitted 
to I?7T° only. 
Data and fits for 
XV, with each 
panel showing 
the indicated 
value of W. Data 
and curve shapes 
and colors are as 
in Fig. [2] The 
curves in all 
panels use the 
same amplitude 
parameters from 
TableQ] 
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Table 1: Results of the fit using two 7 = Breit-Wigner line shapes fitted to only the 2 n channel. The uncertainties 
reflect the stability of repeated fits with varying initial values. "N/A" means no free parameter was allowed. 



Amplitude 


Centroid m R 


Width r 


Phase AO/ 


Flatte y 




(MeV) 


(MeV) 


(radians) 


Factor 


7 = 


1329 + 10 


20 ± 10 


N/A 


1.5+0.3 


7 = 


1390+ 10 


174 + 20 


-0.2 ±0.3 


N/A 



u 300 
| 200 




1 = Amplitude (narrow) 



1 = Amplitude (broad) 



2.2 2.4 2.6 2.8 

W (GeV) „.,..,., 



Figure 4: (Color online) For two 7 = am- 
plitudes fitted to lPn° only, the figure shows 
the strength of each of the isospin amplitudes, 
Co{W), as a function of W, per Eq.[5] 



Next we consider the best fit obtained when using two 7 = amplitudes and a single 7 
I amplitude fitted to the data for all three charge final states: £ + 7r~, iPn and Yrn + . The fit 
parameters are given in Table [5j the reduced^- 2 for this fit was 2.60. Figure [5] shows a sample 
result in one 100 MeV wide bin of W at 2.3 GeV. The colored curves show the total fits to each 
charge state, and the solid, dashed, and dotted black curves show the cross sections that would 
arise from each of the three amplitudes separately. The dotted black curve shows the 7=1 piece 
of the reaction. Its width is similar to the others (54 MeV), with a centroid at 1367 MeV, and a 
Flatte factor that is close to unity, like the one for the 7 = piece. The colored curves show fairly 
successful modeling of the split between the three charge combinations. 



Table 2: Results using two 7 = and one 7=1 Breit-Wigner line shapes and fitting to all final charge states simultane- 
ously. The uncertainties reflect the stability of repeated fits with varying initial values. "N/A" means no free parameter 
was allowed. 



Amplitude 


Centroid m« 


Width r 


Phase AO/ 


Flatte y 




(MeV) 


(MeV) 


(radians) 


Factor 


7 = (low mass) 


1338 + 10 


44+10 


N/A 


0.94 + 0.20 


7 = (high mass) 


1384+ 10 


76 ± 10 


1.8 + 0.4 


N/A 


7 = 1 


1367 + 20 


54+10 


2.2 ± 0.4 


1.19 + 0.20 



Again, the fit was done not just to the W bin shown, but to all nine bins at once. The full data 
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Figure 5: (Color online) Fits of the isospin am- 
plitudes 7 = (solid and dashed black curves) 
and 7=1 (dotted curve) at the given value of W. 
Data and fitted curves for the three charge states 
are: red curve and circles Z + 7r~, blue curve and 
squares: lPn°, green curve and triangles: I,~n + . 
The error bars are combined statistical and sys- 
tematic point-to-point uncertainties. The three 
vertical dotted lines at the left are thresholds for 
the respective decay modes, and vertical dotted 
lines mark the nominal 1 .405 GeV mass and the 
location of the NK threshold m t hi- es h ■ The inco- 
herent background is shown as the thin dashed 
lines. 



set is shown in Fig. [6] In each panel the black curves have the same shapes, differing only in 
overall weight. The way in which the normalizations Cj(W) of the three amplitudes vary with W 
is shown in Fig. [7] Note that it is the square of these coefficients that fixes the relative strength of 
the contributions. The dominant contribution comes from the 7 = amplitude, as expected. The 
7=1 amplitude is weaker in magnitude by a factor of about three, and is about equal in size to 
the smaller of the two 7 = amplitudes. 

It is of interest to compare Table [2] with the previous simpler case shown in Table [T] Upon 
adding the interference-causing 7=1 amplitude, the position and width of the dominant 7 = 
amplitude rises from 1327 MeV to 1338 MeV; both are close to threshold, and agree within 
the uncertainties. The Flatte factor remains in the vicinity of 1.0. The weaker 7 = amplitude 
stays at the same mass, within uncertainties, and gets narrower. Finally, one sees that the 7=1 
amplitude has a centroid at 1367 MeV with a width of 54 MeV. It had a fitted Flatte factor as 
well, and the fit placed it also close to unity. 

Choosing the 7=1 amplitude to have a resonant Breit-Wigner form was, of course, an 
assumption made to create a model. Alternatively one might construct a non-resonant scattering 
amplitude with isovector character. But if this analysis is accurate, and if it can be confirmed that 
an 7 = 1 component in the reaction is reasonably fit with a Breit-Wigner resonance-like structure, 
it raises the question of how to interpret the result. It may point, for example, to the existence of 
a low mass J F - \ state as discussed in Ref. Ir20l . A full theoretical treatment is now needed. 

4. Conclusions 

The model of the mass distribution for the En line shape introduced in Ref. [ 1 1 has been de- 
scribed briefly. It has been applied to recent photoproduction data from CLAS in two alternative 
combinations in order to further explore possible interpretation of the experimental results. 

It has been shown that the 7 = contribution was qualitatively consistent with the picture 
of two poles, modeled as two Breit-Wigner resonances with large channel coupling to the NK 
final state. It has been shown, moreover, that an 7 = 1 component of the reaction mechanism 
can give a reasonable explanation of the variation among charge states in the "En system seen in 
photoproduction. This component was also modeled as a Breit-Wigner line shape with Flatte 
coupling to the NK final state. As shown in Table p] this 7=1 strength is centered near 1367 
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Figure 6: (Color 
online) Fits to 
two / = am- 
plitudes and one 
1=1 amplitude, 
with each panel 
showing the 
indicated value of 
W, without the 
underlying pure 
isospin states 
shown. Data and 
curve shapes and 
colors are as in 
Fig.0 
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MeV. If this is the correct way to model this component of the reaction mechanism, it would lend 
support to existence of a low mass X* state with J p — | . 

Nevertheless, it must be pointed out that a better fit (in the sense of both x 2 and the visual 
appearance of the fit results) was found [ 1 1 previously with an alternative choice of one 7 = 
amplitude and two 7=1 amplitudes. The^ 2 in that study was 2.15. Neither is the model unique, 
nor is the data precise enough to allow only a unique result from the given model. Clearly, more 
modeling and more theoretical work on the interpretation should follow. 
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